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Some examples of the most recent applications of the NRQCD factorization ap- 
proach to quarkonium phenomenology are presented. In the first part of the talk 
the NLO calculations for \cJ an d Tf decays rates are compared to the data and 
the results critically analyzed. In the second part, I show how information on the 
non-perturbative matrix elements can be extracted from the hadronic fixed-target 
experiments and eventually used to test their universality. 



1 Introduction 

A large number^] of achievements, developments, and still-open questions enrich 
the research activity on Non-Relativistic QCD; it would certainly be hopeless 
to cover even a reasonable portion of them here. Fortunately, there are several 
reviews availablea that we can all enjoy and I urge the interested reader to refer 
to them. In this talk, I will therefore concentrate on some specific examples 
which will give us the opportunity to see the factorization approach of NRQCD 
at work. In so doing, it will hopefully be easier to appreciate its strengths and 
to learn how, whenever possible, to avoid its weaknesses. The results presented 
in the following sections are mainly taken from Ref. □ and can be considered 
as applications of the next-to-leading order (NLO) calculations performed in 



2 Inclusive Quarkonium Decays 

The interest for studying quarkonium decays is two-fold. On the one side we 
are faced with a lot of experimental data on masses, total widths and branch- 
ing ratios. On the other, we have at hand predictions for these quantities, in 
most of cases with a NLO accuracy. For instance, the annihilation of quarko- 
nium into final states consisting of leptons, photons and light hadrons can 
provide useful information on the strong coupling at a scale of the order of the 
quarkonium mass M. In general we can write the inclusive decay width as an 

"Invited talk presented at 5th Workshop on QCD, Villefranche-sur-Mer, France, 3-7 Jan 
2000 p. 
b It will suffice to mention that the seminal paper by Bodwin, Braaten and LepageEJ has by 
now almost reached five hundred citations! 
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expansion in the NRQCD operators 

r (Jf = £ f(Q - x){H\0(Q)\H), (1) 
Q 

where T(Q — > x) describes the short-distance decay of a QQ pair in the color, 
spin and angular momentum state Q = 2S+1 L^j' s \ while (H\0(Q)\H) are the 
non-perturbative matrix elements (MEs). Roughly speaking, they represent 
the probability for the heavy-quark pair in the quarkonium H to have quantum 
numbers Q when the annihilation takes place. Even if not calculable within 
perturbative QCD, the MEs have definite scaling properties with respect ta 
the relative velocity in the center-of-mass frame of the two heavy quarks, v.U 
Therefore the sum in Eq. (Q) has to be considered as a simultaneous expansion 
in a s (Mii) and v, and can be truncated at any given order of accuracy. 

There are several sources of uncertainties in these studies. The first one 
is related to non-perturbative MEs that, at least for colour-octet states, are 
poorly known. In fact, the idea of using universality to gain information on the 
MEs is almost unfeasible in inclusive decays. The main reason is that while in 
production one can find processes where the contribution from a given channel 
is enhanced such that values of particular MEs may be obtained (e.g., the 
high-pT J/ip production for (0\Of ( 3 5i)|0)), this is not possible in fully inclusive 
decays. This is the main motivation for studying also semi-inclusive decays, 
where in principle. ihis kind of technique works (e.g., the radiative decays of 
Bottomonium QO'u'lI) . Another source of uncertainty is the renormalization 
scale dependence of the results. As is well known, theoretical predictions at 
fixed order depend on the arbitrary renormalization scale [ir, so for the short- 
distance coefficients we write 

t(Q^x)=A Q a s (» R ) Pa [l+B Q (» R )a s (» R )+0(a 2 s )] (2) 

where pq identifies the leading order, Bq is the NLO correction and Q is 
summed in Eq. ([!]). Unless pg = 0, the coefficient Bq depends both on the 
scale \ir and on the renormalization scheme. The natural expectation is that 
the more terms are added to the perturbative expansion, the less dependence 
on the arbitrary scale remains, and the stability of the result improves. 

In the following sections, I will discuss two cases of interest: the charmo- 
nium P-wave decays and the T hadronic decay. In the first case, we will see 
that, thanks to the measurements available and to the symmetries of NRQCD, 
we can give a satisfactory description of the \cj decays and even formulate a 
prediction of the total width of the h c (the 1 P\ charmonium state). In the sec- 
ond example, we will instead realize that the same strategy applied previously 
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does not work for S'-waves and the data alone are not enough to constrain 
all the unknown parameters. Nonetheless, we will find that if the NRQCD 
expansion in the relative velocity v works for bottomonium (as it should), 
then the color-octet effects in T hadronic decays are of the same order as the 
color-singlet ones and cannot be neglected. 

2.1 Charmonium P-wave inclusive decays 

One of the first applications and successes of the factorization approach of 
NRQCD has been the analysis of the decays of the P-waves.ll3 For such states 
the colour-octet component is at the same order in the v expansion as the 
colour-singlet one: 

\ Xc j) = 0(l)\ 3 P^) + 0(v)\ 3 sf ] g), (3) 

\h c )=0(l)\ipP)+0(v)\is£ ] g), (4) 

and both Fock components have to be considered in computing physical quan- 
tities. The necessity of including all the components at a given order in v to 
obtain a finite result is a general feature of NRQCD. This is due to the fact 
that the non-perturbative matrix elements characterized by quantum numbers 
conserved by the NRQCD interaction at a given order in v mix under renor- 
malization group evolution. To leading-order accuracy in v we can use the 
heavy-quark spin symmetry, which relates the matrix elements: 

{O x )m A c = {xcjPi{*Pj)\xcj) = (hdO^Pi^hc) + 0(v 2 ) , (5) 
(0 8 )ml = (xcAOsCS^lxcj) = (h c \O s CS )\h c ) + 0(v 2 ). (6) 

Moreover, using the vacuum-saturation approximation, which holds true up to 
corrections of 0(v 4 ), we can relate the non-perturbative matrix elements in the 
electromagnetic decay to the hadronic ones. As a consequence, the hadronic 
and electromagnetic decay widths of Xcj, at leading order-in v and at NLO in 
a s , can be read from expressions in Appendix C in Ref . □ and rewritten in a 
compact form as: 

T( Xc j - LH) = Ajna 2 . (l + ^Bj) (O,) + na 2 s (l + ^d) (O s ) , 



(J = 0,2) (7) 

r( Xc i LH) = A 1 a 3 s (O l ) + na 2 s (l + ^D) (0 8 ) , (8) 

L(XcJ - 77) = Ay (l + ^By) ne^aUOi) , ( J = 0, 2) (9) 

L(Xci^77) = 0, (10) 
T(h c -» LH) = A h al(O l ) + na 2 s C h (l + ^-D h ^j (O b ) , (11) 
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Process 



r(xco 



r(xc 2 
r(xc 2 
r(xc 2 
r(xc 2 
r(xc2 
r(xc2 



r(xco 
r(xco 
r(xco 



77) 
77) 



77) 
77) 
77) 
77) 
77) 
77) 



LH) 
LH) 
LH) 



r(xci - LH) 



T( Xc2 -> LH) 



P (McV) 



(4.0 ±2.8) x 10~ 3 
(1.7 ±0.8) x 10- 3 



(0.321 ±0.095) x 10" 3 
(0.7 ±0.3) x 10~ 3 
(3.4 ± 1.9) x 10~ 3 
(0.44 ±0.08) x 10~ 3 
(1.02 ±0.44) x 10~ 3 
(1.76 ±0.62) x 10- 3 



13.5 ±5.4 
15.0 ±4.9 
14.3 ±3.6 



0.64 ±0.10 



1.71 ±0.21 



Reference 



CBALLE 1 
CLEO0 



E76oEi 
CLEO0 
TPC20 
E8350 
L30 
OPAL0 



CBALLtt" 
E8350 
BES 







E760li 



E760E3 



Table 1: Most recent experimental results on Xcj decay widths. The 
errors have been obtained by combining in quadrature the statistical 
and systematic errors given in the quoted references. The widths to 
light hadrons have been obtained from total widths by removing the 
contributions of known radiative decayp-i(Br(xci — * fJ/ip) = 27.3 ± 
1.6% , Br( X c2 -> tJ/VO = 13.5 ± 1.1%). E3 



where A, B, C, D are known numerical coefficients which in general also depen- 
dent onI = 2b Q log and we used the short-hand notation a s = a^ f5 '(/Xfl). 
We note that the hadronic width of the Xd is mainly due to the colour-octet 
component. The result of a global fit of the experimental data (Table [l]) 
to the theoretical expressions for the \cj yields the set of free parameters 
a s , (Oi), (Os), summarized in Table|2[ We verified that the extracted value of 
a s , evolved using the two-loop (3- function, is compatible with the LEP average 
value (a s (Mz) = 0.119). E2l From Table ^ we can read off our best value for 
the non-perturbative singlet matrix element]^ (m c = 1.5 GcV): 

( X c./|Oi( 3 Pj)|Xcj)/m2 = (3.2 ± 0.4) ■ 10~ 2 GeV 3 . (12) 



c For the sake of comparison (with other results, the singlet matrix element is given in the 
original normalization of Ref.tJ. (<D\ = -^j-) 
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m c 


1.5m c 


2m c 


a s (10%) 


0.320 


0.293 


0.273 


{Ox) (10%) 


2.67 


2.41 


2.30 


(Os) (20%) 


2.35 


1.93 


1.84 


d.o.f. 


15.2/10 


15.0/10 


16.1/10 



Table 2: Results of the fit to the non-perturbative MEs and to a s 
for different values of the renormalization scale. Values in MeV. The 
statistical errors from the fit are indicated in parenthesis. 



The above value may be compared with the one obtained from the Buchmuller- 
Tye potential^! (x c j|d( 3 Pj)|Xcj) /m 2 c = 4.8 ■ 10~ 2 GeV 3 and can be useful 
in determining the cross sections for \cj production. As shown in Table || the 
agreement between theory and data improves (i.e., the % 2 /d.o.f. is smaller) 
after including the NLO corrections and, within the assumed approximations, 
the NRQCD approach appears to work quite well even for a heavy-quark as 
'light' as charm. 

Sic stantibus rebus, we can push the analysis further and use the informa- 
tion gathered from \ c j decays to obtain a prediction for the hadronic width 
of the elusive 1 P 1 charmonium state, the h c . In order to curb the errors, we 
used the results to a second fit on the XcJ decay widths with a s fixed to the 
LEP average value. Using Eqs. (||) and Eq. (|ll]), we obtain the prediction for 
the hadronic decay rate of the h c : 

T(h c — > LH) = (0.72 ± 0.32) MeV , (13) 

where the errors from the fit, from the unknown-higher order corrections in v, 
and from the renormalization-scale dependence have been combined in quadra- 
ture. A prediction for the total width of the h c can also be given if the radiative 
width is included: 

Ttot = T(^ c -» LH) + T(h c -» 1Vc ) . (14) 
We can estimate the radiative width by assuming spin symmetry and using §1 

/ E h \ 3 

r( 1 P 1 ^ 7 1 5o)=(^J TfPj^^Sx). (15) 

Plugging numbers!^ into the above expression one indeed verifies an approxi- 
mate independence of J and finds T(h c — > jr/c) ~ .45 MeV. The result for the 



5 



total decay rate 

T TO T(h c ) = (1.2 ±0.3) MeV (16) 

is somewhat higher, but .still compatible within errors with the experimental 
upper limit given in Ref.E3, i.e., r^ c OT < 1.1 MeV (90%CL). We look forward 
to the new data on P- waves that williippefully be collected in the near future 
by the E835 experiment at Fermilab.O 

2.2 T inclusive hadronic decay 

In the case of S'-waves, the spin-singlet states rj and the spin-triplet (vector) 
states J/ip, ip' and the Y(nS), the structure of the Fock state at order v 4 is 

It?) = 0(l)\QQ[ 4 S l o ] }) + 0(v)\QQ[ l Pf ] ]) 

+ 0(v 2 )\QQ[ 3 sf ] }) + O^lQZffsn ) + o(v 6 ) , (17) 

2 

\V) = 0{\)\QQ[ z S [ ?]) + 0{v)\QQ[ 3 P l f ] ]) 

+ O^JIQSpsf]) + 0(v 2 )\QQ[ 3 S?' s] ] } + 0(v 6 ) . (18) 

As we see, in contrast to the P-waves, colour-octet states contribute only at 
order i> 4 , and, at least for bottomonium, where v 4 ~ 0.01, we might expect 
that these effects should be safely neglected. This conclusion is certainly true 
for spin-singlet states when the magnitude of the short-distance coefficients is 
considered: the LO contributions for 1 Sq , 1 S , g 8 ' and 1 p] 8 ' are all at a 2 s . There 
is no enhancement of the short-distance coefficient and so color-octet effects 
for rjcfi can be safely neglected. Different is the case of spin-triplet states: in 
fact, the LO contribution to the hadronic decay of a vector colour-singlet state 
is at order a, while the octets decay already at order a 2 s . We will then include 
these terms, check their relevance and eventually verify whether constrains on 
colour-octet matrix elements can be obtained. 
The total decay width can be written as: 

T(V -> LH) = c i eS^ ] )(V\0 1 ( 3 S 1 )\V) + d^S^iVfP^SJlV) 

+ c 8 ( 3 sf ] )(V\O s ( 3 S 1 )\V) +c s ( 1 S [ * ] )(V\O 8 eS )\V) 

+ J2cs( 3 P l f ] ){V\O s CPj)\V) . (19) 
j 

Using the symmetries of NRQCD one can reduce the number of indepen- 
dent MEs from seven to four. Nevertheless, it is clear that too many MEs 
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are present at order v 4 to be able to extract detailed information from the 
experimental decay rate. An additional constraint can be obtained from the 
width into two leptons, to which only the color-singlet state contributes: 

T(V -> l+r) = c H ( 3 sf)<y|0i( 3 -Si)|F) + du^S^iVlVi^S^V) , (20) 



where cu is now known even at NNLO.BE3 The conclusion is that using the 
available data, we can at most obtain a constraint on a linear combination of 
three color-octet matrix elements. The relativistic correcting to the singlet 



have been calculated numerically inc3 and analytically inc3a and we make 
use of their results. Using the equations of motionEJ to relate the operator O 
to the V, and using heavy-quark spin symmetry to reduce to one the unknown 
MEs of the octet P-waves, we obtain: 

r»<v^-^-.,(.-^l^) 

+ ^<0 8 ( 3 Si)} + ^ ((Os? So)) + 7<0 8 ( 3 P„))) ,(21) 

where the short-hand notation (0({S,P})) = {V\0{{S, P})\V) /m^ has 
been introduced. Some comments are in order: 

• As anticipated, the ratio of the singlet short-distance coefficient to that 
of the octets is typically Ci/cg < cx s /ir and the enhancement is quite 
effective. 



The relativistic corrections are very large for charmonium (~ 100%) and 
relevant for bottomonium (~ 25%). 



At NLO, the corrections for the octets are listed in Appendix_C of Ref.Q. For 
the singlet, the short-distance coefficient can be found in Ref.El. To study the 
effects of the NLO corrections we present some results for a specific choice of 
the adjustable parameters. We fix Aq CD = 237 MeV and /ua = A*R = frrib 
with / = 1/2, 1,2 (from top to bottom), m b = 4.8 GeV and My = 9.46 GeV. 
The result for the LO expression in Eq. (El]) is the following: 

f 11 ) _ ( 0.32] 

T LO (V —>■ LH) = ^ 5.5 } ■ lO^^O^ 3 ^)) + I 0.20 \ (0 8 ( 3 S'i)) 

((O 8 ( 1 5 ))+7(O 8 ( 3 P )}), (22) 
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while for the NLO expression we get: 



Y NLU {V — > LH) = < 4.0 } • 10- J (Oi( J Si)} + I 0.26 ) (O s ( s Si)) 





(O s eS )) + { 6.6 WC 8 ( 3 Po)> • (23) 



As a result, the short-distance-coefficient's variation with different choices 
of the renormalization scale is significantly reduced at NLO. On the other 
hand, the improvement in our theoretical predictions confirms that the typi- 
cal ratio between color-singlet and color-octet short-distance coefficient is ex- 
actly of the same order of magnitude as the expected 0(v A ) suppression be- 
tween the non-perturbativc matrix elements. In other words, unless unnatural 
cancellations between the various octet contributions take place,^ the color 
octet contributions are of the same order as the color singlet one! This quite 
unexpected result leads to the conclusion that the available extractions of 
a s (Q 2 ~ Mr) from the measured ratios of decay widths E3 ( such as = 
r(T hadrons)/r(T -> and i? 7 = r(T -» 7 + X)/T(T -> fj, + ^-) ), 

which neglect color-octet effects, should be regarded with great caution. (See, 
e.g., the discussion on color octet effects in the photon spectrum in the decay 
T — » 7 + X in Ref. □) . This is certainly one of the cases where it would be 
most useful to have good estimates of the non-perturbative matrix elements 
from lattice calculations. 



3 Charmonium Production and Universality 

In complete analogy with Eq. (Q) , we can write the cross section for producing 
a quarkonium state H as a sum of terms, each of which factors into a short- 
distance coefficient and a long-distance matrix element: 

da(H + X) = Y, da(QQ[Q] + x) (O h (Q)) . (24) 

a 

The matrix elements (O h (Q)) are related to the probability for a QQ 
state with definite quantum numbers to evolve into the quarkonium H . In 
contrast to the decays, the large variety of processes in quarkonium production 
at both fixed-target and collider experiments offers the concrete possibility to 

^Remember that at NLO both short- and long-distance coefficients become scheme depen- 
dent and so, in general, non-perturbative matrix elements might not be positive definite. 
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ME 



J/ip i>' Xco 



(0 1 ( 3 Si)> 1.16/6 0.76/6 

(of{ 3 P ))/ml • 10 2 4.0/6 
(0f ( 3 5i)) ■ 10 2 1.06 0.46 0.32 

Table 3: Independent MEs for charmonium production used as input 
values in the fit. Other MEs can be obtained by spin-symmetry relations. 
Values in GeV 3 . 



thoroughly test the factorization approach and in particular the universality 
of the non-perturbative MEs. In fact, NRQCD teaches us that we can trade 
(up to 0(v 4 ) corrections) the color-singlet matrix elements extracted from 
the decays for the ones which enter in the production. Unfortunately this 
is not feasible for color-octet matrix elements. Nonetheless, by studying more 
exclusive processes, as the pt spectrum and/or the polarization of the produced 
quarkonia at hadronic colliders, it is possible to constraint single color-octet 
matrix elements in the NRQCD sum of Eq. (p4[). Eventually we are just left 
with one or two unknown parameters which can be constrained quite well from 
the data. 

In this section we briefly compare the results for hadroproduction of char- 
monium at fixed-target experiments with our predictions. Our analysis follows 
the lines of the one presented in Ref. c3: in particular we use the same experi- 
mental data and just investigate the effects of the inclusion of NLO corrections 
on the extraction of the MEs. _ 

In the following we have used the low-Q 2 MRS A set of parton densities ^ 
rn c = 1.5 GeV ; the renormalization, factorization and NRQCD scales set 
equal to \iu = /if = = 2m c ; the strong coupling, a s , tuned to the one used 
by the PDF set, i.e., A| CD = 152 MeV; and the MEs collected in Table | as 
input values in the fits. Within the above choices the only free parameter left 
is the linear combination of MEs 

M% = (0|Of C5 )|0) + fc (O|0 ^ (3 2 Pj)|O) , (25) 

with k = 7 for the LO analysis and k ~ 6.4 at NLO. The results of the fit are 
shown in Fig. |l| and give: 

,mx>,i a m 3 D [8K _ / 1.8- 10- 2 GeV 3 (J/VO 



M 6A ( b , Fj ) - i Q 2g _ 1Q _ 2 ^3 ^ (26) 
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Figure 1: Left: Total ip' (left) and ip (right) production in proton- 
nucleon collisions (xf > only). For the ip, the total cross section 
includes the radiative feed-down from \ c j and ip' . The solid lines are 

obtained with M^ >A = 2.6 • 10" 2 GeV 3 and M$ A = 1.8 • 10" 2 GeV 3 for 
ip' and ip respectively. 



to be compared with the LO results of E3 

M Lo f i,[8] 3 P [8K _ f 3.0 • 10- 2 GeV 3 (J/V») f27 , 
C*o ,iV o.52 ■ 10- 2 GeV 3 ty'). ( j 

As expected, the inclusion of NLO corrections significantly reduces the values 
needed for the color-octet MEs. Unfortunately, the same sources of uncertain- 
ties present at LO, which are both theoretical (the strong dependence of the 
results on the choice of the heavy-quark mass and on the scale (which is only 
slightly improved at these energies)) and experimental (nuclear dependence, 
limited Xf range, presence of elastically-produced charmonia) make it difficult 
to associate reliable errors to the fitted values. 

Nevertheless, since most of these uncertainties should cancel in the ratios 
of cross-sections, we have also compared 

RX ~^2) (28) 

where 

a tot (J/¥) =a dir (J/^) + Br(^' - JfrP) <r{&) + ]T Br( Xc ./ - Jft) a{ Xc j) , (29) 

J 

with the available experimental data (see Ref. for details). These are shown 
in Fig. |[ While the color-singlet approximation for is rather uncertain and 
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Figure 2: Ratio of direct to total J/ip (left) and Xci/Xc2 (right) produc- 
tion in proton-nucleon collisions as a function of the beam energy. The 
experimental value at 300 GeV for <r dir ( J/ip)/a tot ( J/0) is 0.62 ± 0.04. 



clearly disfavored by the data, both the LO and the NLO NRQCD predictions 
show a much better agreement. Moreover, the NLO result turns out to be 
quite stable under scale variations. With a certain confidence, this result can 
be considered a remarkable success of the NRQCD approach. 

The above analysis shows that, if on the one hand the NRQCD approach 
gives a rather consistent description of the data, on the other hand it is rather 
difficult to extract precise information on the color-octet matrix elements from 
the measurements. In this respect, it is interesting to compare some of the 
extractions available in the literature with the one above. The relevant color- 
octet matrix elements for the ip and ip' are summarized in Table |J Even if 
the very basic assumptions of the cited works are often quite different (PDF 
choice, inclusion of initial-state radiation and/or kx effects, implementation of 
Altarelli-Parisi evolution, etc.) and considering the theoretical uncertainties 
previously mentioned, the overall picture is reasonably consistent. 

4 Conclusions 

I have presented some applications of the NRQCD formalism to the phe- 
nomenology of inclusive quarkonium decays and production. In T decays 
we have established the importance of color-octet contributions to the total 
hadronic decay widths. We have then argued that the present extractions of 
a s from T decays have to be regarded with great care, at least until mea- 
surements of the relevant non-perturbative matrix elements on the lattice are 
available. In the other two examples, the % c j decays and the charmonium 
production at fixed-target experiments, we have been able to extract infor- 
mation on the non-perturbative MEs directly from the data. In the former 
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Process 




ef. 


(^( 3 ^i)) 


k 


Mt 


<Of( 3 Si)> 


Mt 


Tevatron 






0.66 


3 


6.6 


0.46 


1.8 








1.06 


3.5 


4.4 


0.44 


1.8 








0.3 


3 


1.4 


0.14 


0.33 








0.4 


3.5 


6.6 


0.36 


0.78 








1.2 


3.5 


4.5 


0.5 


1.9 


Fixed-target 


EE 


1- 


(0.66) 


7 


3.0 


(0.46) 


0.52 


hadropr. 


* 




(1.06) 


6.4 


1.9 


(0.44) 


0.28 


B- decay 


41 




(1.06) 


3.1 


1.5 


(0.44) 


0.6 



Table 4: Color octet matrix elements for i/> and ip' production as ex- 
tracted from various analysis (* = this work). Values in 10 -2 GeV 3 . 
Numbers in parenthesis are taken as input values. 



case, this has allowed us to make a prediction for the hadronic width of the 
h c charmonium state while in the latter case, we checked the universality of 
the MEs against the available results. The conclusion is that, considering the 
'light' mass of the charm, the NRQCD approach works surprisingly well and 
the extracted MEs are of the order of magnitude predicted by the scaling rules. 
On the other hand, the uncertainties present both at the theoretical and ex- 
perimental level prevent us from performing more precise tests of the most 
distinctive prediction of NRQCD, i.e., the universality of the non-perturbative 
matrix elements. To overcome the above limitations, more exclusive quantities 
are under study. The most exciting and promising one is the expectation that 
high-py. q.uarkonia produced at hadronic-colliders should be transversely po- 
larized. EJC3 Preliminary data from CDFc3 do not confirm this prediction and 
may indicate tkat quarkonium physics is ready for new, unexpected surprises 
(see, e.g., Rcf. o and references therein). 
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